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Abstract 

Semi-inclusive deep inelastic scattering off the Deuteron with production of a slow nucleon in 
recoil kinematics is studied in the virtual nucleon approximation, in which the final state inter- 
action (FSI) is calculated within generalized eikonal approximation. The cross section is derived 
in a factorized approach, with a factor describing the virtual photon interaction with the off-shell 
nucleon and a distorted spectral function accounting for the final-state interactions. One of the 
main goals of the study is to understand how much the general features of the diffractive high 
energy soft rescattering accounts for the observed features of FSI in deep inelastic scattering(DIS). 

Comparison with the Jefferson Lab data shows good agreement in the covered range of kine- 
matics. Most importantly, our calculation correctly reproduces the rise of the FSI in the forward 
direction of the slow nucleon production angle. 

By fitting our calculation to the data we extracted the W and dependences of the total cross 
section and slope factor of the interaction of DIS products, X, off the spectator nucleon. This 
analysis shows the XN scattering cross section rising with W and decreasing with an increase of 
Q^. Finally, our analysis points at a largely suppressed off-shell part of the rescattering amplitude. 
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I. INTRODUCTION 



In recent years, a process that has garnered a fair amount of attention is the d{e, e'ps)X 
reaction at high Q^. In this reaction, deep inelastic scattering (DIS) occurs on a constituent 
of the deuteron and a slow spectator proton is detected in coincidence with the scattered 
electron. This reaction can be used in several ways to study the role of the QCD dynamics 
at nucleonic length scales. At very small spectator proton momenta, the DIS occurs on a 
nearly on-shell neutron and it allows one to extract information about the "free" neutron 
structure function F2N in a way that minimizes the nuclear effects inherent to a bound 
neutron. Detailed information about the neutron structure function helps to constrain the 
QCD models of the nucleon and can be used to determine the relative dtou quark densities 
at large Bjorken x. At larger spectator momenta, high density configurations of the deuteron 
will occur in which the proton and neutron are in very close proximity to each other. Under 
these circumstances the partonic structure of nucleons could strongly modify JT\ with the 
possibility of two nucleons merging into six quark configurations at asymptotically large 
relative momenta in the deuteron [21 E]. Consequently, experiments that explore these 
kinematics can be used to study the modifications of nucleon properties and the role of quark 
degrees of freedom in these situations. Two recent Jefferson Lab Hall B experiments have 
studied the d{e,e'ps)X reactions: one at high [4J, and the other at low spectator momenta 
[S]. New measurements will be possible after the 12 GeV upgrade of JLab is completed. 

In experiments exploring the partonic structure of the nucleon, one generally wants to 
have kinematics that minimize the final-state interactions (FSI) of the produced X-states 
with the spectator nucleon as this FSI make the extraction of the observable one is looking 
for less straightforward. On the other hand the d{e, e'ps)X reaction in kinematics that favor 
larger contributions from FSI can be used in order to study the process of hadronization. 
The attenuation of the produced hadronic state by the spectator when compared to the 
free process can yield information on the space-time structure of the hadronization process. 
Thus in this respect FSI becomes very important part of the semi-inclusive DIS process. 
To quantify the effects of FSI in DIS, model calculations are needed and this has already 
resulted in the development of several theoretical approaches [H [614T2] . 

The major problem one faces in calculations of FSI of DIS products with the spectator 
nucleon in d{e, e'N)X reactions is the lack of the detailed understanding of the composition 
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and space-time evolution of the hadronic system produced after the deep inelastic scattering 
of the virtual photon off the bound nucleon. Moreover both the composition and space-time 
evolution are function of the Bjorken x and probed in the reaction. 

In this paper we study the question on how much the final state interaction of the DIS 
products are defined by the general properties of soft reinteractions. In other words, how 
far we can go with the description of FSI without knowing the specific properties of the 
hadronic intermediate state after the initial DIS scattering? Based on the general properties 
of the reaction a factorized approach is used in the calculations, whereby the cross section 
is split into the parts describing the interaction of the virtual photon with a bound nucleon 
and the distorted spectral function which includes the effect of final-state interactions. The 
deep inelastic interaction with moving bound nucleon is calculated within the virtual nu- 
cleon approximation while the FSI are included using the framework of generalized eikonal 
approximation (GEA) [T314T5] . 

The paper is organized as follows. In Sec. [Tijwe describe the general properties of the reac- 
tion and main assumptions based on which we derive the plane-wave impulse approximation 
and final-state interaction parts of the scattering. An overview of the various approxima- 
tions used in this derivation is also given. In Sec. |III[ the results of our model calculations 
are discussed and compared to the data from the Deeps experiment performed at JLab [1]. 



Finally, conclusions are given in Sec. IV 



II. FORMALISM 



A. General Structure of the Reaction 



We consider the process 

e + d-^e +ps + X, (1) 

in which incoming electron e has energy E^, while Eg/ and dg denote the energy and scattering 
angle of the final electron e'. We define the lab frame four-momenta of the involved particles 
as Pd = {Mo, 0) for the deuteron, q = [i>,q) for the virtual photon (with the z-axis chosen 
along q), Ps = {Eg = ^p^ + m'^,Ps) for the spectator proton a.ndpx = {Ex,Px) = {i^ + Md — 
Es,q — ps) the center of mass momentum of the undetected produced hadronic system X. 
We can express the differential cross section for process ([T| through the four independent 
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DIS structure functions in the following form: 



O"^ 9 \ 

-^ + tan2^cos0F^^(x,g2) + cos20F^2.(x,Q2) . (2) 

Here, c^em is the fine-structure constant, —Q^ = u"^ — (f is the four-momentum transfer, 
Bjorken x = ^ (with m„ the mass of the neutron), y = and is the angle between 
the scattering (e,g) and reaction {q,Ps) planes . 

We now define the nuclear electromagnetic tensor as 

^D=y^\Yl E {DsD\J^''\Xs,,psSs){Xs,,p,Ss\r\DsD) 



47rMz)3 ^ 



X (27r)^5^(g + po - Ps - Px)d-'T^ , (3) 



with (i^r^ a phase-space factor for X, and sd, Sg, and s^; the spin projection of the deuteron, 
spectator proton and X respectively. The four deuteron semi-inclusive structure functions 
FflxyQ"^) are related to components of the nuclear electromagnetic tensor W^'^ as follows: 



FEix,Q') = u^^W^'{x,Q' 
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FEL{x,Q')cos<j) = -2u^W'D%x,Q'' 

m 

F^^ix, Q') cos 20 = u^^iWS'ix, Q') - W^^x, Q')) . (4) 
B. Main Approximations 

In the further derivations we use the following approximations which are based largely on 
the general properties of DIS scattering as well as properties of the subsequent small angle 
rescattering of the fast moving hadronic system off the slow recoil nucleon: 

- virtual nucleon approximation: To treat the electromagnetic interaction with the 
bound nucleon in the deuteron we use the virtual nucleon approximation (VNA) in 
which it is assumed that the virtual photon interacts with the off-shell nucleon in the 
deuteron while the second nucleon is on its mass shell [U [TJ [16] . The VNA is based 



on the following main assumptions: (i) only the pn component of the deuteron wave 
function is considered in the reaction, (ii) the negative energy projection of the virtual 
nucleon propagator gives negligible contribution to the scattering amplitude, and (iii) 
interactions of the virtual photon with exchanged mesons is neglected. Assumptions 
(i) and (ii) can be satisfied when the momentum of the spectator proton is limited to 
Ps < 700 MeV/c [17], while (iii) is satisfied at large (> i GeV^) US HI]. 

The electromagnetic tensor of the interaction is off-shell and the gauge invariance 
is restored by expressing the longitudinal component of the electromagnetic current 
through its O'th component as follows: 

J3 = ^J°. (5) 

The nuclear wave function in the VNA is normalized to account for the baryon number 
conservation [T91422] : 

j a\^D{p)?d^P = l. (6) 

where the light cone momentum fraction of the deuteron carried by 

the bound nucleon normalized in such a way that the half of the deuteron momentum 
fraction corresponds to a = 1. Because of the virtuality of interacting nucleon it is 
impossible to satisfy the momentum sum rule at the same time. As a result 

j a'\^M?d'p<l, (7) 

which can be qualitatively interpreted as part of the deuteron momentum fraction 
being distributed to non-nucleonic degrees of freedom which are unaccounted for within 
the VNA. 

By applying the VNA for calculation of the matrix element {Xsx-,PsSs\J^\Dsd)i we 
can limit the Feynman diagrams taken into account to those of Fig. [l] in which [TJ^a) 
represents the plane-wave impulse approximation (PWIA) diagram. Here DIS occurs 
on the neutron and the proton is left in the on-shell positive energy state without 
further interaction in the final state. The diagram of Fig. [l|^b) shows again DIS on 
the neutron, which is afterwards followed by a X'p — )■ Xp rescattering. In calculating 
this diagram we have to sum over the all possible intermediate X' states. 

The calculation of the final-state interactions is based on the following main assump- 
tions for the rescattering diagram of Fig. [l|^b). 




FIG. 1: Diagrams entering in the model for the d{e,e'ps)X reaction. Panel (a) shows the plane- 
wave contribution. Panel (b) shows the FSI term. 

- diffractive form of the rescattering amplitude: In the considered reaction the FSI rep- 
resents a small angle rescattering of the DIS products off the slow spectator nucleon. 
It is in principle a very complex problem to account for the details of the interaction 
of the intermediate "X"' state since its structure depends on the i^Bj) and the 
produced mass W of the 7*A^ reaction. However in the limit where the produced 
intermediate and final masses are small compared to the transferred momenta: 

g>Mx',Mx', (8) 

one can assume that the propagation of the produced hadronic system is eikonal and 
the general structure of the small angle rescattering is diffractive. The approximation 
of Eq.(|8| allows one to model the FSI amplitude of the hadronic X' system in the 
following form: 

, ^ 2 2 2 B(Q'^,Xgj) 

2_^fx'N,XN = fxN{t,Q ,XBj) =crtot{Q ,XBj){i + e{Q ,XBj))e 2 , (9) 
X' 

where the sum of the all possible X'N — )■ XA^ amplitudes are represented in the 
effective diffractive amplitude form, fxN{t,Q'^,XBj) with effective total cross section 
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(Ttot, real part, e and slope factor B. A similar approximation is used for the FSI studies 
in semi-inclusive DIS scattering [H EHU] as well as for studies of color transparency 
phenomena in which the intermediate state represents an off-shell coherent composite 
system with reduced interaction cross section (see e.g. [131 123VI28] ). In principle, a 
more elaborate model which sums the contribution of different resonances as e.g. in 
Ref. [29] could be used but this would go beyond the goal in this paper of describing 
the reaction with the basic elements of high-energy rescattering. 

factorization: In the situation in which momentum transfer in DIS exceeds the mo- 
mentum of the recoil slow nucleon one can factorize DIS scattering from the amplitude 
of the final state interaction. Such an approximation commonly referred as distorted 
wave impulse approximation (DWIA) is valid in the limit of \/Q^ ^ Ps in which 
case the electromagnetic current is insensitive to the momentum of the stuck nu- 
cleon. The validity of the DWIA was checked quantitatively for quasielastic scattering 
in the case of d{e,e'N)N reactions [T71 130] . These calculations demonstrated that 
for = 2 — A GeV^ factorization approximation works reasonably well for up to 
Ps = 400 MeV/c and then at larger momenta it systematically underestimates the 
FSI contribution as compared to the prediction based on an unfactorized calculation. 
The underestimation can be understood qualitatively, since in the case of nonfactor- 
ization the amplitude of electromagnetic interaction enters in the FSI amplitude at 
smaller values of bound nucleon momenta and therefore predicts more rescattering 
than the DWIA does. This pattern one also expects to be generally valid for inelastic 
interactions. 

approximate conservation law of high energy small angle scatterings: In the eikonal 
regime of small angle scattering there is an approximate conservation law for the "— " 
component [2Z] of slow nucleon momenta involved in the scattering [TBJ. According to 
this law, because the fast particle attains its momentum after the small-angle scattering 
the slow nucleon will conserve its "— " component. This follows from the conservation 
of the "-" component of the total momentum in X'N' — )■ XN scattering and relations 
Px'- ~ "^^'tq^'^ ^ -'■ ^^'^ Px- ^ ^^Vq^^ ^ provided that the condition of Eq. (8) 
is satisfied. This yields: 



Ps'- - Ps- = Px- - Px'- ~ . (10) 



Using this relation and assuming that 

pL < kl (11) 
where is the average transferred momentum in the rescattering one obtains: 

"^i- = {Px' + Ps' - Psf ~ m\, - 2px'±(Ps'± - Psi.) -k]_^ nix, + k]_ > m\, (12) 



where in the above derivation we used the fact that in the hmit of Eq. (11) px'± = 
~Ps'± ~ The above result qualitatively means that in the situation in which two 
collinear particles are produced by the diffractive scattering of a fast and slow particle 
with equal and opposite transverse momenta the mass of the final fast particle is larger 
than the initial mass. 

Using the characteristic values of the diffractive slope, B = 4 — 6 GeV~^, one can 



estimate k±^RMS ~ 500 — 600 MeV/c. This estimate of k± and Eq. (11) further 
constrains the values of spectator nucleon momenta for which the calculations will be 
valid. 

Our derivations in the following two subsections are based on the above assumptions. 



C. Plane- wave impulse approximation 

Applying Feynman diagram rules (see e.g. Ref. [15]) and introducing the effective wave 
functions of the final hadronic system X, the amplitude of the PWIA diagram in Fig. [T]^a) 
takes the following form: 

{Xs^,PsSs\r\DsDf^''' = -^X(PX, Sx)T^^*x-2 o ■ HPs, Ss)rDNN " • (13) 

^ Pi - 

Here, is a wave function for X and F^^.j^- represents the electromagnetic vertex of the DIS. 
The transition of the deuteron into a pn system is described by the vertex function Td^n 
and denotes the spin wave function of the deuteron. The lab frame four-momentum of 
the struck neutron pi in the PWIA is defined as 

p^ = {Md- Es,-Ps). (14) 

We now split the initial nucleon propagator in on-shell and off-shell parts by adding and 
subtracting an on-shell energy part: 

p^ + mn = + + {Ef - ET)l' , (15) 
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with Ef = Md- Es and E""" = ^ml+pl Next we write 

+ m„ = ^M(pi,Si)u(pi,Si) , (16) 

Si 

pofT iron 

{Ef - EDl' ^ \^ ^ 7° E "(P- ' (17) 

where in the last equation we used Si)u{pi, Si) ~ 2mnl, which is consistent with 

neglecting the negative energy component of the bound nucleon propagator. Now, with the 
definition [HI |32] 



we can write Eq. ( 13 ) as 

(fpoS iron \ 
1 + ^- ^7° 
ZrUn J 

x^'}^{piSi,PsS,). (19) 



Even though Eq. (13) is gauge invariant, one can not calculate the off-shell part of the 



current explicitly since the form of the electromagnetic vertex ^t!!*^ x unknown. Instead, 



7* 

in the reminder of the derivation the term — associated with the off-shell behavior 



of the photon-neutron interaction will be dropped and the gauge invariance will be restored 
through Eq. (|5|. Inserting Eq. (19) in Eq. ^ we obtain for the PWIA contribution 



X {2nYS\q + pi-p,)£T,i2nf2Es^^'°ipiSi,Psms)<il'^{pis[,Psm^ . (20) 
We can simplify this further by using 

E ^^'^(PiSi,p,m,)^"«(piS^,p,m,) = E 1"^"° {PiSi,PsSs) fSs,,s[ (21) 

and 

E \'^'°{P^S^ = +l,PsSs)\^ = E l^'^'fe^i = -'^^PsSs)f ■ (22) 
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Eq. (20) then becomes 



After defining the spectral function 



(23) 



(24) 



and using the following expression for the nuclear tensor of the DIS process on a moving 
nucleon 



W 



N 



1 1 



47rm„ 2 



Pi - Px)d^Tx ■ (25) 



X 



we can write Eq. ( 23 ) as 



(26) 



where we also used ~ 1. Substituting Eq. (26) in (4) allows us to relate the four 
deuteron DIS structure functions to the nucleon structure functions. After straightforward 
calculations one obtains the following relations for the deuteron structure functions: 



[Oii + 



+ cos5) -F2N{,Cii,x,Q ) sin 5Fij^iai,x,Q ) 



X S{p,){27if2E, 



F?{x,Q^) = f2Fi;v(«^,£,g') + ^i^27v(«i,£,g')') S{ps){2Txf2E,, 



2, p\ sin 5 



V mi 



F2N{a,,x,Q^)S{p,){2n)'2Es 



Pt 



F^^{x,Q') = 2{1 + cos 6) ^{ai + 



V 



F2N{a„x,Q'')S{p,){2nf2E, 



(27) 
(28) 

(29) 

(30) 



where a,- = a„ = ^ri-, v = — , x 



^,cos(5 = siv?5 



gp, and FiN,F2N are 



the effective nucleon structure functions, which are defined at x and in principle could be 
modified due to the nuclear binding (see e.g. Ref. P). 

Note that the inclusive F2, Fi and structure functions of the deuteron can be obtained 
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from the above given semi-inclusive structure functions through the following relations: 



2,D 



D 



E 

N 

E 

N 



d^Ps 



(27r)22E, 



N 



d^Ps 



(27r)22E, 



d^ps 
2 (27r)22E, 



Fl" = FoD - 2xFi, 



D 



E 

N 



fE + {S^—-x)fP 



'2|gp 



d^Ps 



J (27r)22E, 



d^Ps 



(27r)22E, 



(31) 



where one sums by the contributions of both the proton and neutron. The L.H.S. parts of 
the equations represent the expressions in the case of the Bjorken limit with x fixed and 



D. Final-state interaction amplitude 



With the same notations as in Sec. |II C[ we can write for the amplitude of the FSI 
diagram in Fig. [l|b) 

FSI _ sr^ f d^Ps' ^x{px,sx)u{ps,Ss)Fx'n,xn[P^, + rup] 



GiPx')T';,^,[p., + mn\TDNN ■ X'^ 



(32) 



[Px' - ^x' + ^e] [pI -ml + ie] 
where G{px') describes the Green's function of the intermediate state X' which has four- 
vector px' = Pi' + q = (i^ + Md — Esi,q — Ps') and a mass m^', while the intermediate 
struck neutron has four- vector pii = {Mc — Es',—ps')- Fx'n,xn represents the invariant 
X'N — )■ XN scattering amplitude which is expressed in the following form 

Fx'N,XN{s,t) = a/(s - (m„ - mx'y){s - (m„ + mx'Y)fx'N,XN{s,t) 

= (3{s, mx')fx'N,XN{s, t), , (33) 

with s = {px + Ps)"^ = {px' + Ps'Y the total invariant energy of the scattering system 
and the scattering amplitude fx'N,XN defined such that lm[fx'N,XN(t = 0)] = atot, where 
(Jtot represents the total cross section of the scattering of the produced X' system off the 



spectator nucleon. Based on the assumptions of the VNA from Sec. II B, the intermediate 
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spectator nucleoli can be placed on the nucleon mass-shell by integrating dPps' through the 
positive energy pole only: 

. (34) 



This allows us to use the on-shell spinor relation + nip = ^u{ps',Ss')u{ps',Ss') in 



the nominator of Eq. (32). For the propagator of the initial neutron we again use the 



prescription of Eqs. (15) to (16). For the intermediate state X' an on-shell relation for the 



Green's function G{px') = J2s ,i'{Px',Sx')4'''{Px',Sx') is used as in the high limit the 



off-shell contribution in Eq. (32) becomes small due to the large momentum involved in the 
propagator of the intermediate state X'. By making use of Px = {q + Pd — PsY = ^x, the 
denominator of the X' propagator can be rewritten as 

Px' - ^x' + ie = 2\q \ {ps',z - Ps,z + A + ze) , (35) 

with 

A = I ^. (Es - Es.) + ^ I ^.^ . 36 



All this combined with the deuteron wave function of Eq. (18) allows us to write the FSI 
amplitude as 

{Xsx,PsSs\r\DsDf''' = -Y. E 77^,Pis,mx'){Xsx,PsSs\fx 'N,XN{s,t)\X'Sx',Ps'Ss') 

X' Si,s^,,s^, ^ ^ 

AE,.\q\[p,.^z-p,^z + A + ie] VI ; . I J 

In a next step, we assume the rescattering amplitude conserves the heli cities of all particles 
involved 

{XSx,PsSs\fx'N,XN{s,t)\X'Sx',Ps'Ss') ~ {X S^.PsS s\f X' N,Xn{,S ,t)\X' S^.Ps' S s)5s„s^,5s^,s^, , 

(38) 

and we use the following approximation to take the current matrix element out of the 
integration: 

{X'Sx\T^,j^^X'\Pl'^i) ~ {XSxlV^^N^xlPi^i) ■ (39) 



This allows us to factorize the nuclear tensor again like in Eq. (26). For the sum of the 
plane-wave and FSI amplitudes, we then obtain 



W^^" = W^''S'''^'-{ps){2nf2E, , (40) 
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with the distorted spectral function defined as 



\dist. 



1 ^'^{Pis,,psSs) -Y, j 



(27r)3 4 I g-l 



x{XS:j:,PsSs\fx'N,XN{s,t)\X'Sx,Ps'Ss) 



'^''''{pi'Si,Ps'Ss 



[Ps',z - Ps,z + A + ie] 



(41) 



E. Distorted spectral function 



For calculation of the distorted spectral function in Eq.(41 ) we use VNA model of deuteron 



wave function of Eq. (18) which can be related to the non-relativistic deuteron wave function 

(42) 



^d{p) 



D 



D 



2{Md - E,) ' 

which explicitly conserves the baryonic sum rule of Eq. (|6]). The parameterizations for 

the non-relativistic wave function used in this paper all take the following form (see e.g. 
Refs. [331 EH): 



*i°(pSl,-pS2) 



E 



p^ + m" 



p2 _|_ 



X 



sd 



(43) 



with S{p) 



1 ( 3(71 ■pa2 -P 



p2 ^ — c?i ■ (T2 j the tensor operator. Such form allows to perform the 



dps'^z integration in the distorted spectral function of Eq. (41) analytically by making use 



of the pole structure of these parameterizations as well as the pole of the propagator in 
Eq. (41) at pgi^z = Ps,z — A. In the latter case the mass of the produced intermediate state 



mx' enters in the phase factor A. We note that even though we sum over the all possible 
intermediate states X' the mass rrix' is defined by the four momenta of the interacting 
virtual nucleon and virtual photon q. Based on the assumption that FSI is dominated by 
a small angle diffractive scattering, the phase factor A is evaluated based on the property 
of the approximate conservation law of "— " components of rescattering particle momenta 



discussed in Sec, II B Taking into account the relation of Eq. (12) in the definition of the A 



factor in Eq. (36) we evaluate: 



A 



{Es - rup) + 



rrix-l 

2|g| 



A = I ^1 — [Es - nip) 



for 7 < m\ , 
for 7 > m\ , 



(44) 
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where 7 = m\,{pii = 0) = + 2m„z/ — Q"^ is the produced DIS mass off the stationary 
nucleon. The latter approximation for m^' is justified by the fact that due to the peaking 



of the deuteron wave function at small momenta the integrand in Eq. (41) is dominated by 
smaller virtual nucleon momenta than in the PWIA term 



The Psi,z integration in Eq. (41) is performed analytically by closing the integration 



contour into either the upper or lower complex hemispheres. In both cases [38] one obtains: 



dps', 



' Ps',z - Ps,z + A + ie 



-iTT'^'"'{Ps,z,Ps',±,Si,Ss) - TTPs,z'^{Ps,z,Ps',±,Si,Ss) (45) 



where the distorted wave function \E' is defined in Eq (47). In the above equation the first 



term can be identified with the (imaginary) on-shell part of the ps'^z propagator while the 



second term with the (real) principal value integration. Inserting Eq. (45) into Eq. (41) one 
obtains for the distorted spectral function: 



S{p, 



\dist. 



SD,Ss,Si X' 



d'^Ps',± (3{s,mx') 
(27r)2 4 1^-1 ^MEs 



X [{X,p,\f^%^^^is,t)\X',p,)^^-ip,,, 



-i{X,Ps\f^N,XN{s,'t)\X',Ps')Ps',z'^"'iPs',Si,Ss) 



where for the distorted wave function of the deuteron one obtains 

'ttu I T) _ I 'tnrx I r» . I r'^ 

^'^{Ps,Si,S2 



Mps) + — T^^yPs) + — 7^ — T- {<s{ps) - S{ps,±)) 



PIz 



(46) 



(47) 



with 



Mp) = Yl 

j 

wi{p) = 



p\ + m|p2 _|_ 772^ 
dj 



J ■^z pj_ + m'jp'^ + mj 



W2 



(p) = E 



dj 



sjpi + m) 



(48) 



The distorted spectral function in Eq.(46) depends on the intermediate state X' through 



the final state rescattering amplitude only. As a result one can factorize the sum over the 
X' in the form of Eq.(|9]) and represent the on-shell forward scattering amplitude in the form 
of 



rxN = (^UQ\W){i + e{Q\W))e 



(49) 
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with W the invariant mass of the produced hadronic state X. For the off-shell amplitude 
there is no clear prescription, but following our main goal of studying the semi-inclusive 
DIS based only on basic properties of the high-energy scattering we identify two extreme 
cases for off-shell part of the rescattering amplitude, one when it is taken to be zero {no off- 
shell FST) and the other in which off-shell amplitude is assumed to be equal to the on-shell 
amplitude f^j^ referred as maximal off-shell FSI. The numerical estimates for we used 
in our calculations will be discussed below in Sec. IIIII 

With this, we have all the ingredients needed to compute the cross section of Eq. ([2j). 



III. RESULTS 



A. Experimental observables 



In this section, we compare calculations in our model with the first results extracted from 
data taken in the Deeps experiment at JLab [3]. Events of the data set were binned in Q^, 
Ps, cos 6s (with Og = q^s) and x (or the invariant mass of the produced hadronic state W). 
In order to compare our model calculations with the data, we integrate Eq. ^ over </)e', use 



dx 
dx 



2x^ V 



1 

2£ 



X \q 

and relate FiAr(aj, x, Q^) to F2Ar(ai, x, Q"^) for a moving nucleon: 
FiAr(ai,x,Q^ 



2x 



1 + R 



+ 



an 



2x 



Pt 



2Q2 



R 



F2iv(ai,x,(5^) 



(50) 



(51) 



0.18 is the ratio of the longitudinal to transverse cross sections for scattering 



where R 

off the nucleon. This yields for the differential cross section: 
da 47rQ;EM 1^1 A x^y^^i\ I 2tan^ 



dxdQ'^d^Ps 



Q'^x rUr, 



1-y 




l + R 



aj 1 
— + 



2x 



-1 



+ 



2Q2 



F2N{ai,x,Q^)S{ps). (52) 



Now using Eq. ( [52| ) we need to reproduce the quantity F2nP{Ps) (with P{ps) = 
2{Md-^Es) I^S^(P'^)I^) which the experimental data are given in Ref. j^. For this we 



divide the cross section of Eq. (52) with the following prefactor 

pix^y"^ 1 



Ana 



EM 



Q^x 



y 

2(1 + i?) 



+ (1 - y) + 



l + R, 



(53) 
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where y = This results to the following representations of our model calculations: 

(F2ArP)model = ^ f j-j^2^S ) " (^^) 

In numerical estimates we use the SLAC parameterizations for the neutron structure 
functions Fi^ and F2N [35] in the calculations as these were used in the analysis of the 
Deeps data [1]. The arguments of the nucleon structure functions are defined from the off- 
shell kinematics {pi + q) = W^, where the four momentum of the initial nucleon is defined as 
Pi = Pd — Ps- No additional modifications due to nuclear modifications like the EMC effect 
are assumed for the nucleon structure functions. This is in accordance to our approach of 
estimating the properties of the reaction based on the basic properties of the high energy 
scattering rather than modeling the specific details of the reaction. Additionally we deem 
the influence of these modifications small in comparison with the typical magnitude of the 
experimental uncertainties to extract unambiguous information here. 



B. Numerical Estimates 



We start first with the calculation of the quantity of Eq. (54) for the typical kinematic 
setting of the experiment ^ with = i g GeV^ VV^ = 2 GeV and Ps = 390 MeV/c. 
For on-shell part of the rescattering XN — )■ XN amplitude we use the diffractive form of 
the parameterization of Eq. ^ with characteristic values of atot = 50 mb, B = Q GeV~2 
and e = —0.5. Our estimate of the total XN cross section is based on the assumption 
that the final state consists of the hadronic state equivalent to one nucleon and one pion. 
In principle it is possible to develop specific model (see e.g. Ref. [12]) describing the XN 
rescattering, however we follow here our main goal of understanding how far we can go with 
describing data on basic properties of high energy scattering. For the off-shell part of the 
XN rescattering we use two limiting cases as discussed above: fxN~^^^^^ — ('^^ off-shell 
FSI) and fxN~'^^'^^^ ~ /^V'*'"^" {maximal off-shell FST). The results of these calculations are 
given in Fig. |2} 

As the figure shows, FSI effects continuously grow in the forward angles of production of 
recoil proton. This result is strikingly different from the case of the quasielastic d{e, e'N)X 
scattering in which case the FSI is maximal at transverse angles (~ 70°) of recoil nucleon 
production (see e.g. [T7j) and diminishes in the forward direction. The continuously growing 
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FIG. 2: (Color online) Comparison between the Deeps data [3| and model calculations at 
= 1.8 GeV^. The dashed black curve is a plane- wave calculation, the other include final-state 
interactions. The effective total cross section and slope parameter in the final-state interaction 
amplitude are fixed to dtot = 50m6, B = 6GeV^^ and e = —0.5. The dotted blue curve has an 
off-shell rescattering amplitude amplitude equal to the on-shell one [maximal off-shell FSI), the 
dash-dotted red curve has no off-shell FSI. 



FSI contribution in the forward direction for DIS scattering follows from the specific structure 



of the phase factors (A) entering in Eq. (44) which follows from Eq. (12). For forward 
angles the dominant mass contribution m'x decreases, as can be seen in the stationary 
approximation (virtual photon energy u decreases with forward angles). As m'x decreases 



the off-diagonal mass term in Eq. ( 44 ) grows bigger and so does A causing the peak to shift 
to more forward angles. 

Another observation from Fig. [2] is the relatively small contribution due to the off-shell 
part of the XN rescattering amplitude. This result is in agreement of the space-time analysis 
of high energy small angle rescattering of Ref. [36], according to which the longitudinal 
distances that off-shell particle propagates before rescattering significantly shrinks in the 
high and fixed Bjorken x limit. This results in the suppression of the off-shell part of 
the FSI. 
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The next question we address is whether the parameters of X'N rescattering amphtude 
are sensitive to the produced DIS mass W and Q^. For this we assume some of the parameters 
entering the rescattering amphtude of Eq. ^ to be free parameters. We made fits using 
one (effective total cross section atot) or two ((Xtot and slope factor B) free parameters. The 
real to imaginary part ratio of the amplitude was fixed at e = —0.5, a value extrapolated 
from nucleon-nucleon scattering parameterizations. 

For the off-shell rescattering amplitude in addition to above mentioned no off-shell FSI 
and maximal off-shell FSI options we consider the third approach in which case we param- 
eterize the off-shell amplitude as 

ffN = Ct(Q^ + 6-(g^ py^))e- ""'t-""' ^ , (55) 

were the effective cross section and real part parameters were taken equal to the on-shell 
ones, but the slope parameter was taken as a new free parameter in the fit. This will give 
us some measure of the size of the suppression as compared to the on-shell amplitude, this 
approach is referred to as fitted off-shell FSI. 

The parameters were fitted for each (Q^, W) to all measured spectator momenta. When 
comparing the results of the fits to the data it became clear that the model fits systematically 
underestimate the data at the highest measured spectator momentum ps = 560 MeV. This 
may be a consequence of the factorization used in this model, which begins to break down 



at these momenta (see discussion in Sec. JIB). In subsequent fits we decided to exclude the 



highest spectator momentum as we deem the model not adequate enough to describe the 
data in these kinematics. 

Figures [3] to [6] show the results of these fits in which both the effective cross section 
a and slope factor B were free parameters. The plane-wave calculations generally show 
little dependence on the spectator angle, in clear disagreement with the change seen in the 
data. The calculations including FSI manage fairly well to describe the data over the covered 
kinematics. When comparing the three off-shell descriptions, we see that differences between 
the three become smaller with higher spectator momentum. This indicates the diminished 
importance of the off-shell part of the rescattering amplitude in these kinematics. 

At the lowest missing momentum of Ps = 300 MeV, there is an oscillating structure in the 
data which disappears for high W but is still present in the calculations. When comparing 
the three calculations including FSI, we see that there is a large difference between them in 
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FIG. 3: (Color online) Comparison between the Deeps data [3] and model calculations at = 
1.8 GeV^ at measured values of invariant mass W and spectator momenta p (= ps in the text) of 
300, 340 and 390 MeV. The dashed black curve is a plane- wave calculation, the other include final- 
state interactions. The effective total cross section and slope parameter in the final-state interaction 
amplitude are fitted parameters for each W, the real part is fixed at e = —0.5. The dot-dashed 
green curve only considers on-shell rescattering, the dotted blue curve has an off-shell rescattering 
amplitude equal to the on-shell one and the full red curve uses the off-shell parameterization of 



Eq. (55). 
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FIG. 4: (Color online) Comparison between the Deeps data [3] and model calculations at = 
1.8 GeV^ at measured values of invariant mass W and spectator momenta p (= Ps in the text) of 
460 and 560 MeV. Graphs as in Fig. [3} 



the backward angles. There, the no off-shell FSI calculation is smaller than the plane-wave 
calculations while the maximal off-shell FSI calculation becomes significantly bigger. The 
fitted off-shell FSI calculations sits somewhere in between and tends to agree more with the 
maximal off-shell FSI calculation at low W and with the no off-shell one at high W. At this 
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FIG. 5: (Color online) Comparison between the Deeps data [3] and model calculations at = 
2.8 GeV^ at measured values of invariant mass W and spectator momenta p (= Ps in the text) of 
300, 340 and 390 MeV. Graphs as in Fig. [3| 



value of spectator momentum, the plane-wave and final-state interaction amplitudes are of 
comparable magnitudes |39]. This makes the final result quite sensitive to small variations 
in the FSI amplitude and its off-shell description, thus providing some way of explaining the 
larger discrepancy between data and different calculations as compared to higher ps values. 
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FIG. 6: (Color online) Comparison between the Deeps data [3] and model calculations at = 
2.8 GeV^ at measured values of invariant mass W and spectator momenta p (= Ps in the text) of 
300, 340 and 390 MeV. Graphs as in Fig. [3| 



At higher spectator momenta, the no off-shell FSI calculations more or less exhibit three 
regimes. At backward angles they almost coincide with the plain- wave calculations. Around 
90 degrees they show a steep rise, which flattens out at the forward angles. This agrees 
with the intuitive picture of final-state rescattering. The maximal off-shell FSI calculations 
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on the other hand have a more constant slope for the whole of the spectator momentum 
range. The calculations with a fitted off-shell FSI description show the best agreement with 
the data, which is to be expected as they have an extra free parameter. Over the whole 
of the kinematics they generally agree more with the no off-shell FSI calculations than 
the maximal off-shell FSI ones, pointing at a largely suppressed off-shell amplitude. At 
the highest measured spectator momentum the FSI curves systematically underestimate the 
data, pointing in the direction of a breakdown of the factorization used in this model. 




FIG. 7: (Color online) The fitted values of effective cross section a and slope factor B for the no 
off-shell FSI calculations used in Figs. [3] to [6] as a function of the invariant mass W. Full blue 
curve is for = 1.8 GeV^, the dashed green curve for = 2.8 GeV^. 

The final question we addressed in the above described fitting procedure is whether the 
data indicate on and W dependence of the parameters of XA^ rescattering. 

Figures [7] and |8] show the values of the fitted parameters Utot and B used in respectively 
the no off-shell and the maximal off-shell FSI calculations. At PF ~ 1.2 GeV (corresponding 
to the production of a A), we get a cxtot around 40 mb. For the higher invariant masses, 
the cross section drops to around 20-25 mb and rises with increasing W, consistent with 
the production of more hadronic constituents in the intermediate state of the DIS reaction. 
The cross section doesn't flatten out at the highest W, showing that hadronization occurs 
before the rescattering in these kinematics. With increasing Q^, the value of the (XN) cross 
section parameter also becomes consistently smaller in this region, indicating reduced final- 
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FIG. 8: (Color online) The fitted values of effective cross section a and slope factor B for the 
maximal off-shell FSI calculations used in Figs. [3] to [6] as a function of the invariant mass W. Full 
blue curve is for = 1.8 GeV^, the dashed green curve for = 2.8 GeV^. 

state interactions. This could be a sign of an onset of a color transparency effect, in which 
with increasing the hadronic state is produced in a state with smaller transverse size, 
subject to reduced QCD interactions with the medium. The values for the slope parameter 
B are also largely correlated with those of Utot with a smaller slope parameter at higher Q'^ 
and larger B for higher W, although we also see some clear deviations from this picture (e.g. 
a.t W = 2.4 GeV in the no off-shell FSI fit). Overall our fitting procedure indicates that the 
availability of more and W data points may allow to gain important insight about the 
and W dependence of the total cross section of NX scattering. 

IV. CONCLUSION 

Based on the virtual nucleon approximation framework, we developed a model to describe 
semi-inclusive deep inelastic scattering of the deuteron. To describe the final-state interac- 
tion of the spectator nucleon with the produced hadronic state X, the general features of 
diffractive soft rescattering were used, without specifying the structure or space-time evo- 
lution of X. The generalized eikonal approximation was used to calculate the scattering 
amplitudes based on effective Feynman diagram rules. A factorized approach was used to 
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split the cross section into a part describing the virtual photon interaction with the off-shell 
neutron and a distorted spectral function containing the final-state interactions. 

The model calculations were compared to data taken in the Deeps experiment at Jefferson 
Lab. We first compared our calculation with the data for typical kinematics of the experi- 
ment with characteristic parameters for final state interactions. This comparison indicates 
a good agreement with the data most importantly describing correctly the rise of FSI in 
forward direction. This result is opposite to what observed in quasi-elastic kinematics. 

To gain insight on the and W evolution of the FSI further calculations were done in 
which two free parameters (effective cross section a and slope factor B) in the rescattering 
amplitude and three different off-shell rescattering prescriptions were considered. Results 
were fitted for each (Q^, W) to the available data. The fitted off-shell rescattering parame- 
terizations yielded results similar to the calculations with only an on-shell rescattering am- 
plitude included over a wide range of the kinematics, giving evidence for a largely suppressed 
off-shell rescattering. The resulting calculations showed reasonable agreement between the 
data and the calculations including final-state interactions. There were some discrepancies 
at the highest spectator momentum, which may be caused by the breakdown of the factor- 
ization used in the model. At the lowest Ps = 300 MeV there is also an oscillating structure 
in the calculations which isn't exactly present in the data at higher W. The calculations in 
this case proved to be very sensitive to the size of the off-shell amplitude. 

When inspecting the values of the parameter fits, three features emerge: i) The effective 
cross section rises with increasing W, consistent with the creation of more hadronic con- 
stituents taking part in the rescattering. ii) There is no evidence for a plateau at the highest 
measured W values, indicating that the hadronic state has hadronized before rescattering 
takes place, iii) We obtain lower values for utot for the higher Q'^ value, which could be 
interpreted as a sign of emerging color transparency. However, more data at higher Q'^ are 
needed to make more definitive statements. 
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